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We report the first observation of the associated production of a W boson and a Z boson. This 
result is based on 1.1 fb _1 of integrated luminosity from pp collisions at ^/s = 1.96 TeV collected 
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interaction are produced individually in large numbers 
in pp collisions at the Fermilab Tevatron, and their cross 
sections have been measured with high precision [lj . The 
production of heavy vector boson pairs (WW, WZ, and 
ZZ) is far less common and can involve the triple gauge 
couplings (TGCs) between the bosons themselves via an 
intermediate virtual boson. Deviations of measured di- 
boson production properties from standard model (SM) 
predictions could arise from new interactions or loop ef- 
fects due to new particles at energy scales not directly ac- 
cessible to a given experiment 0|. At the Tevatron, TGCs 
are probed at the highest energy scales yet achieved. 

In this Letter, we report the first observation of WZ 
production. The production is observed in pp collisions 
at sfs = 1.96 TeV using 1.1 fb _1 of integrated lumi- 
nosity collected by the CDF II detector at the Fermilab 
Tevatron. We consider the decay channel WZ — > £'vg>££, 
where I 1 and £ are electrons or muons directly from W 
and Z decay, respectively, or from the leptonic decay of 
r's when one or both vector bosons decay to r leptons. 

The most sensitive previous search for WZ production 
was reported by the D0 Collaboration using 0.3 ftT 1 of 
integrated luminosity, where three WZ — > i'vull candi- 
date events were found 0]. The observed events had a 
probability of 3.5% to be due to background fluctuations, 
corresponding to a(WZ) < 13.3 pb at 95% C.L. A search 
for the sum of WZ and ZZ production in decays to 2, 
3, and 4 lepton channels by the CDF Collaboration us- 
ing 0.194 fb^ 1 of integrated luminosity determined that 
a{WZ + ZZ) < 15.2 pb at 95% C.L. 0. The next-to- 
leading order (NLO) WZ cross section prediction for pp 
collisions at Js = 1.96 TeV is 3.7 ± 0.3 pb 0. 

The components of the CDF II detector relevant to 
this analysis are described briefly here; a more complete 
description can be found elsewhere @. The detector ge- 
ometry is described using the azimuthal angle (/) and the 
pseudorapidity r\ = — ln[tan(#/2)], where 9 is the polar 
angle with respect to the proton beam axis (positive z- 
axis). The pseudorapidity of a particle originating from 
the center of the detector is referred to as ?7d. 

The trajectories of charged particles (tracks) are re- 
constructed using silicon microstrip detectors @, H| and 
a 96-layer open-cell drift chamber (COT) [§] inside a 
1.4 T solenoid. The number of COT layers traversed 
by a particle in the range |?7d| < 1 is 96 and decreases 
to zero for \r\&\ 2. The silicon system provides cov- 
erage with 6 (7) layers with radii between 2.4 cm and 
28 cm for \r/ d \ < 1.0 (1.0 < \r)d\ < 2.0). Outside of the 
solenoid are electromagnetic (EM) and hadronic (HAD) 
sampling calorimeters, segmented in a projective tower 
geometry, and constructed of layers of lead or iron ab- 
sorber, respectively, and scintillator. The EM section is 
the first 19-21 radiation lengths (^o), corresponding to 
one hadronic interaction length (A) and contains electro- 
magnetic showers, while the HAD section extends to 4.5- 
7 A and contains the majority of a hadronic shower. The 



calorimeters are divided into central (|^d| < 1-1) and for- 
ward (1.1 < |r?d| < 3.64) regions. Outside of the central 
calorimeters are muon detectors consisting of scintillators 
and drift chambers. 

Including the leptonic r decays, the branching fraction 
of the WZ state to three eor/j leptons is 1.8%. When 
coupled with the small SM cross section, this implies that 
only a small number (~70) of WZ — » i'vi'll events are 
expected to be produced in 1.1 fb _1 at the Tevatron. Fur- 
thermore, in order to identify a W Z event in this decay 
channel, all three charged leptons must be detected. The 
CDF II detector, however, has gaps in calorimeter cover- 
age and limited forward (|r7d| ~ 1) tracking efficiency. 

In order to maximize the total acceptance, while min- 
imizing the backgrounds from jets and photons misiden- 
tified as leptons, we exploit all available reconstructed 
tracks and clusters of energy in the EM calorimeter. We 
separate these into seven non-overlapping lepton cate- 
gories: three each of electrons and muons, and a sev- 
enth for tracks that are non-fiducial to the calorimeters 
and are not identified as muons. In this context, "non- 
fiducial" refers to detector regions that are inactive for 
energy measurement because they are either not covered, 
or are only partially covered, by calorimeter components. 

All lepton candidates are required to be isolated such 
that the sum of the Et for the calorimeter towers in a 
cone of AR = y/(Aij) 2 + (A</>) 2 < 0.4 around the lepton 
is less than 10% of the Et for electrons or pt for muons 
and track lepton candidates. The transverse energy Et of 
an energy cluster or calorimeter tower is ii'sin^, where E 
is the associated energy. Similarly, pt is the component 
of track momentum transverse to the beam line. 

All electron candidates are required to have a cluster 
of energy in the calorimeter with the ratio of deposition 
in the HAD to EM sections consistent with being due to 
an electron. These candidates are divided into three cat- 
egories: those in the central calorimeter, those in the for- 
ward calorimeter matched to a track, and those in the for- 
ward calorimeter without a matched track. The central 
electron category requires a well-measured COT track. 
Since the tracking efficiency is low in the large |?yd| region, 
a track pattern algorithm which starts with calorimeter 
information and attempts to attach silicon hits is used 
for forward electrons. For forward electrons without a 
matched track, both charge hypotheses are considered 
when forming WZ candidates, since the charge is deter- 
mined from the track curvature. 

For all muon candidates, the energy deposition in both 
the EM and HAD calorimeter sections is required to be 
consistent with that of a minimum ionizing particle. The 
muon candidates are divided into a category in which the 
tracks match to reconstructed track segments ("stubs") 
in the muon chambers and two categories of tracks that 
do not match to stubs ("stubless"). The stubless muon 
candidates are designated as central or forward, depend- 
ing on the calorimeter to which the track is projected. 
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The stubbed and central stubless muons have strict re- 
quirements on the number of COT hits and the \ 2 of the 
track fit in order to suppress background muons from K ± 
or 71"^ decays. To increase the track finding efficiency in 
the forward region, we use an algorithm that starts with 
silicon detector hits in addition to one that starts with 
COT hits. The forward stubless muons require at least 
60% of the traversed COT layers to have hits. To sup- 
press the background from cosmic rays and K ± or 7r ± 
decays, we require the point of closest approach of the 
track to the beamline to be consistent with having orig- 
inated from the beam, in addition to using a cosmic ray 
rejection algorithm. 

An additional category consists of tracks that neither 
project to the fiducial regions of the calorimeters nor 
are identified as stubbed muons. The requirements for 
these track-only lepton candidates are the same as for 
central stubless muons, but without the calorimeter re- 
quirements. Due to the lack of calorimeter information, 
electrons and muons cannot be reliably differentiated for 
this category, and are therefore treated as having either 
flavor in the WZ candidate selection. If an electron or 
track-only candidate is consistent with a photon conver- 
sion, as indicated by the presence of an additional nearby 
track with a common vertex, the candidate is rejected. 

To measure the presence of a neutrino, we use missing 
transverse energy $ T — \ Et,% nr,i\, where fiT,i is the 
transverse component of the unit vector pointing from 
the interaction point to calorimeter tower i. The $J T 
calculation is corrected for muons and track-only lepton 
candidates, which do not deposit all of their energy in 
the calorimeter. 

The events we consider must pass one of four online 
trigger selections. The events with central electrons re- 
quire an EM energy cluster with Et > 18 GeV matched 
to a track with pr > 8 GeV/c. The events with forward 
electrons require an EM energy cluster with Et > 20 
GeV and an uncorrected, calorimeter-based measure- 
ment of $ T > 15 GeV. Muon triggers are based on 
stubs from the muon chambers matched to a track with 
Pt > 18 GeV/c. Trigger efficiencies are measured in lep- 
tonic W and Z data samples 

The WZ candidates are selected from events with ex- 
actly three lepton candidates using requirements that 
were optimized with Monte Carlo simulation without ref- 
erence to the data. At least one lepton is required to sat- 
isfy the trigger and have E T > 20 GeV (p T > 20 GeV/c) 
for electrons (muons). We loosen this requirement to 10 
GeV (GeV/c) for the other leptons to increase the WZ 
kinematic acceptance. Aside from WZ production, other 
SM processes that can lead to three high-py leptons in- 
clude dileptons from the Drell-Yan Z/7* process (DY), 
with an additional lepton from a photon conversion (Z7) 
or a misidentified jet (Z+jets) in the event; ZZ produc- 
tion where only three leptons are identified and the un- 
observed lepton results in and a small contribution 



from it — > WbWb, where two charged leptons result from 
the W boson decays and one or more from decay of the b- 
quarks. Except for it, these backgrounds are suppressed 
by requiring $ T > 25 GeV in the event, consistent with 
the unobserved neutrino from the leptonic decay of a W 
boson. We also require the azimuthal angle between the 
ifjrp direction and any identified jet with Et > 15 GeV 
or WZ candidate lepton to be greater than 9° to sup- 
press DY backgrounds in which the observed $ T is due 
to mismeasured leptons and/or jets. 

We require at least one same-flavor, opposite-sign lep- 
ton pair in the event with an invariant mass Mg+£- in 
the range [76, 106] GeV/c 2 , consistent with the decay of 
a Z boson. This range is referred to as the "Z-mass re- 
gion." If there is more than one such pair, the leptons 
with Mi+£- closest to the Z mass [l(| are treated as the 
Z boson decay candidate pair. In order to suppress back- 
grounds from ZZ, we require that no additional track in 
the event with pt > 8 GeV/c, when combined with the 
lepton that is not part of the Z boson decay candidate 
pair, has an invariant mass in the Z-mass region. The 
overall acceptance for WZ — > I'l/pil, using the described 
selection criteria, is 13.4%. 

The acceptances for the WZ, ZZ, Z7, and it processes 
are determined using Monte Carlo calculations followed 
by a GEANT-based simulation [ll| of the CDF II detec- 
tor. The Monte Carlo generator used for WZ, ZZ, and 
it is PYTHIA 12| and for Z7 is the generator described in 



[131 ] . We use the CTEQ5L parton distribution functions 
(PDFs) to model the longitudinal momentum distribu- 
tion of the initial-state partons 14|. An efficiency cor- 
rection, of up to 10% per lepton, is applied to the simula- 
tion based on measurements of the lepton reconstruction 
and identification efficiencies using observed Z — > l + t~ 
events. An additional correction is applied to the Z7 
background estimate based on a measurement of the pho- 
ton conversion veto efficiency in data. The background 
from Z+jets is estimated from a sample of events with 
two identified leptons and a jet that is required to pass 
loose isolation requirements and contain a track or energy 
cluster similar to those required in the lepton identifica- 
tion. The contribution of each event to the total yield 
is scaled by the probability that the jet is identified as 
a lepton. This probability is determined from multijet 
events collected with a set of jet-based triggers. A cor- 
rection is applied for the small real lepton contribution 
using single W and Z boson Monte Carlo simulation. 

Systematic uncertainties associated with the Monte 
Carlo simulation affect the Z7, ZZ, it, and WZ simula- 
tions similarly. The uncertainties from the lepton selec- 
tion and trigger efficiency measurements are propagated 
through the analysis, giving uncertainties of 1.2% — 2.0% 
and 0.4% — 0.9% for the respective efficiencies of the dif- 
ferent signal and background processes. The uncertainty 
due to the $ T resolution modeling is determined from 
comparisons of the data and the Monte Carlo simulation 
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in a sample of dilepton events. For WZ, ZZ, and tt pro- 
duction, where it is an observed particle that produces 
the observed $ T , we determined the uncertainty to be 
1%. The uncertainty due to the $J T modeling for the 
Z7 background is much larger (25%) than for other fi- 
nal states because it depends on the non-Gaussian tails 
of the resolution function. The uncertainties on the ZZ, 
Z"i, and it cross sections are assigned to be 10% 0, 10% 



and 15% [III, 17 1, respectively. For the Z7 back- 
ground contribution, there is an additional uncertainty 
of 20% from the detector material description and con- 
version veto efficiency. The detector acceptance varia- 
tion due to PDF uncertainties is assessed to be 2% using 
the 20 pairs of PDF sets described in [l8| . The system- 
atic uncertainty on the Z+jets background is determined 
from differences in the measured probability that a jet is 
identified as a lepton for jets collected using different jet 
Et trigger thresholds. These variations correspond to 
changing the parton composition of the jets and the rel- 
ative amount of contamination from real leptons. The 
total systematic uncertainty on the Z+jets background 
prediction is 23%. The signal and background estimates 
obtained from simulation have an additional 6% uncer- 
tainty originating from the luminosity measurement [l^ | . 

In addition to the signal region ($ T > 25 and M e + e - 
in the Z-mass region), we define two independent trilep- 
ton control regions, both with $ T < 25 GeV, but differ- 
ent Mi+i- criteria, to validate our background estimates. 
The "Z-mass control region" is defined to have Mg+ e - in 
the Z-mass region and is dominated by Z+jets and Z7 
where the photon is from initial-state radiation. The "Z- 
veto control region" is defined to have M e +£- outside of 
the Z-mass region and a minimum value for Mg+g- of 
40 GeV/c 2 . This region is dominated by Z7 where the 
photon is from final-state radiation. We expect 243.5 ± 
38.8 (250.9 ± 48.3) events in the Z-mass (Z-veto) control 
region and observe 215 (241) events. The results for the 
trilcpton classifications we consider are shown in Table U 
and are in good agreement with the expectations. 

The expectation in the signal region for WZ and each 
background contribution is summarized in Table [TXJ In 
this region, we expect 2.7 ± 0.4 background events and 
observe 16 events. A breakdown of the observed (ex- 
pected) events by flavor classification is as follows: 6 (2.7 
± 0.2) eee, (2.0 ± 0.2) eefi, 1 (1.5 ± 0.1) e/i/z, 1 (1.2 
± 0.1) fififi, 5 (2.0 ± 0.2) eet u 2 (1.3 ± 0.1) efi£ t , 1 (1.1 
± 0.1) fi/dt, (0.5 ± 0.1) el t l u and (0.2 ± 0.1) fj,£ t £ t . 
Here It denotes the track-only lepton candidates having 
unknown flavor. The distributions of $! T , Mg+i- , and the 
W transverse mass M^f = ^2Et$t(^ ~~ cos </>£„), where 
<\>i v is the azimuthal angle between the non-Z candidate 
lepton and the Ifi T direction are shown in Fig. [1] The 
data are in good agreement with the SM prediction. 

Due to the unobserved neutrino, events from W Z — > 
Vvfjti are expected to have larger $ T on average than the 
DY and ZZ backgrounds. We exploit this information 



TABLE I: Summary of the expected and observed yields in 
the trilepton control regions. In the classification column, l t 
denotes a track-only lepton candidate having unknown flavor. 
The eee, e/u/z, eeit, and e\itt classifications receive a large con- 
tribution from Z7 events where the photon is reconstructed 
as a forward electron without a matched track. 



Flavor 


Z-mass 




Z-veto 




Classification 


Expected 


Data 


Expected 


Data 


eee 


116.3 ± 19.2 


103 


114.8 ± 22.5 


103 


e e fi 


1.8 ± 0.3 


2 


1.4 ± 0.4 


4 


e /i fi 


62.5 ± 10.3 


50 


69.2 ± 14.0 


62 


fi n n 


1.1 ± 0.2 


1 


0.3 ± 0.1 


1 


e e It 


29.6 ± 4.6 


20 


33.5 ± 6.2 


31 


e (j, £ t 


24.9 ± 4.1 


33 


26.5 ± 5.2 


34 


H H it 


2.7 ± 0.4 


5 


1.9 ± 0.4 


3 


e l t It 


4.0 ± 0.7 


1 


2.6 ± 0.5 


2 


H It It 


0.4 ± 0.2 





0.4 ± 0.1 


1 


Total 


243.5 ± 38.8 


215 


250.9 ± 48.3 


241 



TABLE II: Expected number of events in the signal region 
for W Z and the background contributions. "Lumi" refers to 
the integrated luminosity uncertainty, which is absent for the 
Z+jets because it is determined from the same dataset. 



Source 


Expectation ± Stat ± Syst ± Lumi 


Z+jets 


1.21 ± 0.27 ± 0.28 ± - 


ZZ 


0.88 ± 0.01 ± 0.09 ± 0.05 


Z 1 


0.44 ± 0.05 ± 0.15 ± 0.03 


tt 


0.12 ± 0.01 ± 0.02 ± 0.01 


Total Background 


2.65 ± 0.28 ± 0.33 ± 0.09 


WZ 


9.75 ± 0.03 ± 0.31 ± 0.59 


Total Expected 


12.41 ± 0.28 ± 0.45 ± 0.67 


Observed 


16 



by performing a binned maximum likelihood fit for the 
signal yield using the following $ T bins: 25 < $ T < 45 
GeV and $ T > 45 GeV. This binning was chosen to 
maximize our sensitivity for a WZ observation without 
reference to the data. We expect 2.0 ± 0.4 (0.7 ± 0.1) 
background events and 6.5 ± 0.5 (5.9 ± 0.4) TVZ events, 
with 9 (7) events observed in the lower (upper) $ T bin. 
We define A In £ as the log of the likelihood ratio between 
this fit and the no signal hypothesis. For our data, we 
find 2Aln£ = 37.8. We interpret this result using 10 10 
background-only Monte Carlo experiments, out of which 
only 11 had a larger value of 2Aln£ (a probability of 
l.lxl0~ 9 ), corresponding to a significance equivalent to 
six standard deviations. 

This result represents the first observation of WZ pro- 
duction. The measured cross section is 

a( P p -> wz) = s.ot{;l ( stat -) ± °- 4 ( s y st -) P b > 

consistent with the SM expectation. 
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FIG. 1: Distributions for WZ candidates of (a) the $ T (b) the dilepton invariant mass for the same-flavor opposite-sign dilepton 
pair closest to the Z mass, and (c) the W transverse mass calculated from the remaining lepton and the Ip T . In (a) and (b), 
the arrows indicate the signal region. 
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